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ABSTRACT 

We present a first-stage study of the effect of using knowledge from electromagnetic (EM) observations in the gravitational wave 
(GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in the low- 
frequency range, 10~ 4 Hz < / < 1Hz. In particular, we examine the extent to which the accuracy of GW parameter estimation improves 
if we use available information from EM data. We do this by investigating whether correlations exist between the GW parameters 
that describe these binaries and whether some of these parameters are also available from EM observations. We used verification 
binaries, which are known as the guaranteed sources for eLlSA and will test the functioning of the instrument. We find that of the 
seven parameters that characterise such a binary, only a few are correlated. The most useful result is the strong correlation between 
amplitude and inclination, which can be used to constrain the parameter uncertainty in amplitude by making use of the constraint of 
inclination from EM measurements. The improvement can be up to a factor of ~ 6.5, but depends on the signal-to-noise ratio of the 
source data. Moreover, we find that this strong correlation depends on the inclination. For mildly face-on binaries (( < 45°), EM data 
on inclination can improve the estimate of the GW amplitude by a significant factor. However, for edge-on binaries (( ~ 90°), the 
inclination can be determined accurately from GW data alone, thus GW data can be used to select systems that will likely be eclipsing 
binaries for EM follow-up. 
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1. Introduction 

The space-based gravitational wave (GW) detector in considera- 
tion by ESA, eLISA, is expected to observe millions of compact 
Galactic binaries ( Nelemans 2009] |Amaro-Seoane et al.||2012| l 
with periods shorter than about a few hours, amongst other as- 
trophysical sources, and resolve several thousand of these bina- 
ries ( |Nissanke et al.||2012| ). About 50 compact binary sources 
have been observed at optical, UV, and X-ray wavelengths (e.g. 
|Roelofs et al.|2010 i. The types of binaries known to us are in- 
teracting systems (AM CVn stars, ultra-compact X-ray binaries, 
and cataclysmic variables) and detached systems (double white 
dwarfs (WDs) and double neutron stars (Nelemans 2009 , 201 1 )). 
The AM CVn stars are binary systems where a WD accretes 
matter from a low-mass, helium-rich (hydrogen-deficient) ob- 
ject (Solheim 2010). Their mass transfer is driven by GW ra- 
diation loss. The known ultra-compact X-ray binaries consist 
of neutron stars that are accretors whose donors are inferred to 
be either helium rich or carbon/oxygen rich ( |in't Za nd 2005). 
Double WDs are predicted to be the most common systems, 
which sometimes tend to be the outcome of many binary evo- 
lutionary paths (Webbink 1984). Of all these known systems, a 
handful lie in the eLISA band and will be individually detected. 
These are known as verification binaries since they are guar- 
anteed sources for the detector. Parameter uncertainties in the 
verification binaries and Galactic binaries in general have been 
studied in the literature extensively by using Fisher information 



matrix (FIM) analyses (e g. |Cutler 1998 ; Ta kahashi & Seto|2002 
Stro eer & Vecchio|2006| ). These studies have been done for var- 



ious configurations of classic LISA ( |LISA Study Team] 1998 ), 
which was designed to have a larger baseline of five million km 
with six laser links interchanging between three stations located 
at the vertices of a triangle that was to house two proof masses 
each ( Vallisneri et al.|2008| ). Instead eLISA will have a baseline 
of one million km with four laser links interchanging between 
the proof masses. The parameter uncertainties depend intricately 
on the observation conditions and the geometry of the detector 
( |Takahashi & Seto|2 002 ) and the same is true for the (possible) 
correlations between the parameters of a Galactic binary. In this 
study, we wish to quantify whether any such correlations exist 
that could be useful in constraining the GW parameter estimates. 
Some of the GW parameters are the same as or related to electro- 
magnetic (EM) parameters, (e.g. inclinationcw = inclinationeM, 
/gw = 2/Porb, etc. ). Thus, we can use an independent (EM) con- 
straint of a (EM/GW) parameter that correlates strongly with an- 
other (GW) parameter to improve the accuracy of the latter. This 
work provides a first step towards developing strategic plans in 
observing those potentially useful EM parameters that can im- 
prove the GW parameter accuracy. In this paper, we present a 
FIM analysis to study whether correlations within GW param- 
eters exist. For the useful correlations that we find, we predict 
quantitatively the improvement in the GW parameter when there 
is prior EM data in the correlated parameter. The paper is struc- 
tured in the following way. In Section 2, we briefly summarise 
the signal models, both the instrumental and foreground noises, 
and our data analysis. We present our results and interpretations 
in Section 3. Finally, we discuss how the results differ from the 
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old LISA detector in Section 4 and present our conclusions in 
Section 5. 



2. Signal modelling and data analysis 

2.1. Gravitational wave signals from a Galactic binary 

We consider three verification sources, AM CVn, 
SDSS J0651+2844 (hereafter J0651), and RX J0806.3+1527 
(HM Cnc), whose physical parameters are summarised in 
Table Q] 

AM CVn and HM Cnc are mass-transferring systems known 
to astronomers as AM CVn binary systems, which were de- 
scribed in the introduction. J0651 ( |Brown et al.||20l"T| ) is an 
eclipsing detached WD binary system that was spectroscopically 
identified in the Sloan Digital Sky Survey (SDSS) catalogue. 
AM CVn and J0651 can be modelled as monochromatic GW 
sources, which means that they are described by seven param- 
eters: dimensionless amplitude iJ[), frequency (/), polarisation 
angle (iff), initial GW phase (</>o), inclination (cos i), ecliptic lati- 
tude (sin j3), and ecliptic longitude (A). HM Cnc can be modelled 
as a mild chirper with an additional eighth parameter, the chirp- 
ing frequency (/). The two polarised gravitational waveforms 
used in the strain for slowly evolving binaries are given by (e.g. 
IKrolak et al.|2004) > 



1 + COS 2 L 

h+(f) = M cos 



(2nft + Tift 1 + «/> ) ; 



h x (t) = &l cos i sin (27r/f + nft 2 + <poj . 



where 



and 
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In these expressions, Ai = im\m.2) i ^ Hm\ + m.2) 1 ^ 5 is the chirp 
mass and d is the distance to the source. The monochromatic 
waveforms are given by setting / = in the expressions above. 

2.2. Detector response to the GW signals 

If a GW signal is present, then the output of a detector will con- 
tain the strain, h(t), and noise, «(?)• Thus, the detector registers 



s(t) = hit; 0) + n(f), 



(5) 



where 6 is a vector characterising the seven (or eight) parameters 
of the binary. Since eLISA will have a motion around the Sun and 
a cartwheeling motion around its centre of mass, a monochro- 
matic signal from a WD binary will be modulated in its ampli- 
tude, frequency, and phase in complicated ways. The resulting 
signal will be spread over a range of frequency bins of the de- 
tector ( Cornis h" & Larson"]|2003 1. Galactic binaries typically ra- 
diate monochromatic signals at low freque ncie j*] and thus the 
response at the detector can be written as ( |Cornish & Larson] 
2003) 

hit) = Ait) cos¥(f), (6) 



where 



ah) = KF + (OMor + (F x it)h x it)y] 



2-1 1/2 



(7) 



1 low / is a function of the detector transfer frequency, /,,/ <k /,. 
The /, is defined according to the detector armlength, L, i.e. f, = 
clilTiL). For eLISA, L = 10 9 m, /, w 5 X lfr 2 Hz 



The functions F + X (f) are the antenna beam patterns of the de- 
tector, and they depend on the source's sky position (A, /3), its 
orientation it//), and the detector configuration. The phase of the 
signal is given by 



«P(f) = 2k ft + 0o + O D (0 + <M0, 



(8) 



where <I> D (f), Op(f) are frequency (Doppler) and phase modula- 
tions, respectively (Cornish & Larson 2003). Doppler modula- 
tion is given by 



D (f) = InfL/c sin/? cosi2nf m t - A), 



(9) 



where f m — 1 /year is the modulation frequency. The phase mod- 
ulation is given by 



<M0 



• arctan 



F + h. 



(10) 



Different architectures of the triangular space-based interferom- 
eter generate a number of independent data streams that provide 



different responses to the incoming GW signal ( Vallisneri et al. 
2008 ). For the most recent interferometer design in considera- 
tion, eLISA, the output is a single unequal-arm Michelson data 
stream, X. This is a linear combination of phase shifts measured 
at the different spacecraft (by comparing the incoming light with 
a local reference source) shifted in time in such a way as to rep- 
resent interference between two light beams travelling through 
the arms of the detector in opposite ways, which is a particu- 
lar implementation of so-called time delay interferometry (TDI, 
|Armstrong et al.| (fl999). Val lisneri] ([2005a) provide a detailed 
description of how TDI works and an explanation of why it pro- 
duces an interferometry signal in which the phase shift induced 
by a passing GW is preserved, while the much larger shifts in- 
duced by instrumental noise are strongly suppressed. We made 
use of the existing numerical software Synthetic LISA (Vallisneri 
2005b) to simulate accurate time-domain series of the instru- 
mental noise and the GW signals in the form of these TDI ob- 
servables. The data stream is a discrete series for a given obser- 
vational time Tob S , where the samples are separated by At. The 
detector response to GWs and instrumental noises have been dis- 
cussed extensively in the literature (e.g.|Cutler|1998[|Cornish & 
|Rubbo||2003j |Kr61ak et al][2004l |Vallisneri||2005bl , hence we 
only summarise the most essential expressions relevant to our 
data analysis. In Synthetic LISA, the strain at the detector, sit), is 
modelled as the TDI X observable. This is the quantity we work 
with in our data analysis. 

2.3. Noise 

There are two types of noise to consider: instrumental noise 
and Galactic foreground noise due to unresolved compact bina- 
ries. For the particular geometry of eLISA used here, the instru- 
mental noises (mostly from shot noise), acceleration noise, and 
other types of noise (e.g optical bench noise), are characterised 
by their power spectral densities (PSD^jt of 2.31 x l0~ 3S f 2 , 
6 x 10~ 48 .r 2 , and 2.76 x 10~ 38 / 2 , respectively, in units of Hz" 1 . 
The instrumental noise is modelled as a random, Gaussian pro- 
cess. We note that the sampling time, At, should be carefully 
chosen while simulating the instrumental noise in order to cor- 
rectly interpret the TDI observables. The choice of At should 
correspond to a frequency that is several times higher than the 



2 Note there is no PSD for laser noise, since we assume that it is 
completely cancelled in the TDI observable X. 
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Table 1. Physical properties of verification binaries based on observations summarised in Nelemans (|201 1\. The signal- 
to-noise ratio (S/N) listed in the second column is from GW data analysis, as explained in Section 3. 
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B.l in Appendix 



highest frequency where the TDI responses have to be analysed 
d Vallisneri||2005b| l . This means that since AM CVn and J0651 
have relatively low frequencies they can be analysed with sam- 
ples of Af = 64s, whereas for HM Cnc, we need a lower sam- 
pling time of at least Af = 16s (see Figure 
B). 

The foreground noise from the Galactic binaries is simu- 
lated using the Lisasolve (Vallisneri 2011) software where ev- 
ery binary (monochromatic and/ or mild chirpeip]l is modelled 
in the frequency domain ( [Cornish & Li ttenberg 2007 ). This dif- 
fers from simulating signals using Synthetic LISA, where the sig- 
nal is modelled accurately in the time-domain. Lisasolve instead 
makes use of the very slowly evolving nature of the binaries to 
approximately model the signals both directly and speedily in 
the frequency domain. We only use double-detached WD bi- 
naries because they form the majority of the foreground noise 
( |Nelemans et al.||200l[ |Petiteau||20T2"l >. We include more than 
2.7 x 10' detached Galactic double WDs from a simulation with 
the same assumptions about binary evolution and Galactic dis- 
tribution as those in Nelema ns et al.| ( |2004[ ), but with (about 
a factor of ten) higher intrinsic resolution. These assumptions 
have also been used for the simulations of the Galactic binaries 
on wh ich the Mock LIS A Data Challenge (MLDC) rounds are 
based ( |Littenberg|20TT) . 

The unresolvable Galactic foreground noise is obtained by 
iteratively subtracting the resolvable sources (i.e. binaries ex- 
pected to be detected) from the simulated population as follows: 



1 . The S/N (see Eq. 14 1 is computed for each binary against the 
Gaussian instrumental noise and initial Galactic foreground 
noise for T D b s of two years. The initial foreground noise is 
calculated using the initial catalogue of 2.7 x 10 7 detached 
Galactic double WDs. 

2. All the sources with S/N > 5 are removed from the initial 
catalogue/dataset. The reduced dataset is used to simulate 
the reduced Galactic foreground. 

3. Using the reduced Galactic foreground and the same instru- 
mental noise, the S/N for each of the binaries in the reduced 
catalogue is calculated. The process is iterated with step 2. 

We applied a perfect subtraction of the bright sources where 
any spurious effects in the data set were not taken into account 
(since the sources were removed from the population before gen- 
erating the signals). A Markov chain Monte Carlo based data 
analysis of the MLDC shows that all the recovered parameters of 
the subtracted sources have a strong peak of zero bias when com- 
pared to their injected values in the training data set (Littenberg 
|201 l| l. Thus, this perfect subtraction scheme is not expected to 
introduce strong biases in our results. From our subtraction pro- 
cedure outlined above, we estimate the number of resolved WD 
binaries to be ~ 11,000. Using the same subtraction procedure 



3 Mild implies that /// <k l/T obs , where T obs is the observational 
time. 



and S/N threshold, Niss anke et al.| ( |20T2) l estimate the number of 
bright sources to be half of our estimate, although their estimate 
is only for a T D b s of one year. Furthermore, using a higher S/N 
threshold of seven, the number of resolvable sources is 3,000 for 
an observation time of two years ( Amaro-Seoane et al. 2012 ). 
If we use a threshold S/N of seven, we find 4,500 resolvable 
sources. The PSD of the unresolvable Galactic background is 
consistent with findings in the literature where the foreground 
noise for eLISA is almost at the level of instrumental noise for 
this detector (Nissanke et al. 2012 ), unlike in the case of the clas- 
sic LISA where the foreground noise was predicted to dom inate 
at / < 3mHz ( |Nelemans et al.|2001| [Timpano et al.|2006 1. 



2.4. Data analysis 

For GW sources with known waveforms, one can use matched 
filtering methods (Finn 1992| Cutler & Flanagan 19941 to ex- 
tract the signal parameters and estimate their uncertainties. 
Consequently, when the noise is Gaussian, the parameter uncer- 
tainties are given by their joint Gaussian probability distribution 
function ( |Cutler|1998| l 

p(o-e) = Vdet(r/27r) exp |- ^ Ty o- 9l cr^j , (11) 

where F is known as the Fisher information matrix (FIM) given 
by 

' dh 



dh 



(12) 



The inner product (...|...) is a generalisation of the time-domain 
correlation product and is conventionally defined as 



(a\b) = 4 df 



a*(f) b(f) 



S n(.f) S n (fo) Jo 



dt a(t) b(t). (13) 



Eq 



1' 



holds for quasi-monochromatic binaries that have an al- 
most constant noise PSD, S„(f), in the frequency region where 
the binary radiates (Cutl er|1998| l. The S/N of a source is defined 
as the inner product of a signal with itself 

S/N 2 = (h\h). (14) 

In the limit of signals with a high signal-to-noise ratio (S/N 
» 1), the inverse of the FIM gives the variance-covariance ma- 
trix C - T l . The diagonal elements C„ give variances (or mean 
square errors) in each parameter, ((cr,) 2 ), and the off-diagonal el- 
ements describe the covariances (or correlations) between them. 
For each of our verification binaries, we calculate this matrix 
to investigate the correlations between the binary parameters. 
The derivative in Eq. 12 dh/dO,, is numerically calculated in the 
time-domain 

h(t ; 0i + dBi) - h(t; - d6 t ) 



h'(t; 0j) = 



2ddi 



(15) 



The latter equality follows from Parseval's theorem. 
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where d9, should be chosen carefully. In general, the quantity 
d6i should be as small as the machine accuracy allows for, but 
not too large to suffer from the truncation erroirl Thus, for well- 
behaved functions, dOi ~ yfe 6 C , where e ~ 10 is the machine 
accuracy and 9 C is some typical value of the corresponding pa- 
rameter ( Press et al. 2002| l. To find a good choice of c/6*,, we com- 
pute <Tj for a range of d9i with logarithmic intervals and select 
the value for each corresponding parameter around which the 
standard deviations become stable. By stable, we mean that in- 
creasing or decreasing d6j by an order of magnitude should lead 
to values of <t, that vary by no greater (smaller) factor than 1 . 1 
(0.9). An example of the stabilisation of the variance-covariance 
matrix is provided in Table [A~T| of Appendix A for the case of 
AM CVn. We perform this stability check for all our verification 
binaries. For instance, in the case of AM CVn, d0j[ ~ 10~ 30 , 
d9f ~ 10~ n , dBt.i/, ~ 10~ 8 etc. All our analysis is done for 
r obs = 2 years. 

3. Results 

Here we list the variance-covariance matrix C for AM CVn, 
HM Cnc, and J0651 and discuss the strongest correlations that 
we found. The off-diagonal elements are specified by the nor- 
malised correlations, Cy, and the diagonal elements are quoted 
as the square root of the variances (or standard deviations), c,-,-, 
i.e. 

C'j i— 
cy = ca = y/Ctt = o~i- (16) 

Thus, Cy can have values in the range [—1,-1-1] where, cy = +1 
means maximally correlated and Cy = -1 means maximally 
anti-correlated. We consider highly correlated parameters 6>,, 8j 
to be those for which |cy| > 0.8. In the later subsections, we 
show that these correlations are affected when the inclination is 
varied. 

The correlation matrices for the three verification binaries 
with observed parameters as shown in Table [T] are listed below. 
The measured GW parameter values are shown above each 
matrix. Since there are no EM measurements for </>o and tf/, we 
set them to n and n/2, respectively, for all three binaries. The 
choice of the </>o value does not influence the results shown 
below, whereas the choice for the if/ has an effect that will be 
addressed later. In the matrices, we take the medians of values 
obtained using 50 different instrumental-noise realisations and 
the standard deviations about these medians (cr(S/N) and <x(<x,)). 
They are listed above and below each of the corresponding 
matrices, respectively. The values of S/N we find here are com- 
parable to the ones found by other groups (Petiteau 2012). The 
normalised correlations are accurate within the quoted precision 
and hence their standard deviations from the instrumental noise 
are not listed. 

J0651,S/N= 10.72 ±0.23. 

Strong correlations (as defined above) are printed in bold face. 
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AM CVn, S/N = 11.54 ±0.19. 
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3 This error comes from higher-order terms in the Taylor-series ex- 
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3.1. The correlation between J{ and cos l 

The amplitude and inclination of AM CVn and HM Cnc (when 
modelled without the chirp, /) are highly anti-correlated, with a 
value of Cjti cosl ss -0.99. The reason for this is that binaries with 
inclinations in the range from 0° to ~ 45° (face-on or close to 
face-on) have signals with the same structure (see Figure [TJ, so 
the signals in this range can be transformed into each other by 
simply changing either the amplitude or the inclination. Thus, if 
we change the amplitude of one of these signals within its GW 
uncertainty, we can obtain a very similar signal by changing its 
inclination within its GW uncertainty as well. However, for in- 
clinations in the range ~ 45° - 90° (either close to edge-on or 
edge-on), this scaling does not apply. This is because the signals 
in this range not only vary in terms of the overall scale but also 
in their structure, when their inclination is changed. This is es- 
pecially evident for a binary with ( ~ 90° (lines of t = 60°, 90° in 
Figure [TJ. Since the inclinations of both AM CVn and HM Cnc 
are within the range where similar signals can be obtained by 
varying either the inclination or the strain amplitude, the param- 
eters cos l and J[ are highly correlated. Furthermore since J0561 
is an edge-on system with ( ~ 90°, both cos l and J[ are distin- 
guishable, as can be seen in Figure[T] This is potentially the most 
useful correlation from the GW analysis where an EM measure- 
ment of the inclination of a binary can constrain the error in the 
strain amplitude M. For example for AM CVn, an EM constraint 
in i with a 1 -<x accuracy of 4° ( Roelofs et al.|2006] > improves the 
uncertainty in J[ from 1 .084 x 10" /z to 1.65 X 10" 23 as shown 
in Figure [2] The corresponding 1-cr accuracy in J[ is estimated 
by selecting the two-dimensional probability distribution func- 
tion (2D PDF) points that lie in the range 41° < i < 45°. This 
improves amplitude accuracy by a factor of 6.5 compared to the 
use of GW data alone. 

HM Cnc's inclination measurement derived from the EM ob- 
servation is quite uncertain compared to that for AM CVn. This 
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Fig. 1. Left: Envelopes of the GW signal (TDI X observable) for binaries with parameters of AM CVn, but for various inclinations 
as a function of time (innermost to outermost envelopes correspond to t = 90°, 60°, 45°, 40°, 20°, 0° respectively). The modulation 
of the signal is due to the complicated annual motion of eLISA around the Sun. For low inclinations (t = 0° - 45°), the signals only 
differ in amplitude, while the higher inclinations also differ in structure. To highlight this, we plot in the right panel the normalised 
ratio of the upper envelopes compared to the l = 0° envelope using the same line styles (but now the t = 90° is the outermost line). 
High inclination systems clearly have a unique structure, while low inclination envelopes only differ in amplitude, illustrating the 
degeneracy between l and 




Fig. 2. Left: Two-dimensional error ellipses of M and cos t extracted from the variance-covariance matrix, C, for AM CVn. The 
two thin black lines represent 1-cr and 2-cr ellipses and the black dashed lines represent the 1-cr GW errors in M and cost with 
the cross at the true value of the parameters. The unphysical values (i.e. negative numbers along the amplitude axis) are caused 
by the Gaussian tails of the parameter uncertainty about their true values. The right panel is a zoom of the area that is constrained 
by the 1-cr EM error of 4° shown in dash-dotted lines. The 1-cr GW standard deviation in J{ decreases from ~ 1.08 x 10~ 22 to 
~ 0.165 x 10~ 22 , roughly a factor of 6.5. The corresponding error ellipses for the reduced PDF are shown as thick black ellipses for 
1-cr and 2-cr. 



is mainly due to the large uncertainty in the mass ratio and the 
assumption about the GW radiation that goes into estimating the 
inclination (Roelofs 'et al.|2010) . Hence, if we take an EM con- 
straint in i with a 1-cr uncertainty of 7°, the amplitude improves 
by a factor of 1.34 as shown in Figure [3] However, if we were to 
obtain a better i constraint of 4° for HM Cnc, J{ would improve 
by a factor of three, which is still a smaller improvement than 
that of the AM CVn's. This is due to the higher S/N of HM Cnc: 
if we place AM CVn at a closer distance such that it has a higher 
SNR of ~ 40, the 1-cr accuracy of 4° constrains the amplitude 



more tightly by a factor of 2.95, which is very similar to the case 
of HM Cnc. By constraining the strain amplitude, we constrain 
the chirp mass, the distance, or a combination of the two ^j- 
(see Eq.[3]l. The GW frequency, the third contributor to the am- 
plitude, is already determined very precisely. Additionally, if / 
is measurable for a chirping binary, then from Eqs.[3]and|4]one 
can constrain the error in the distance ( Schu tz|1996| ). 
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Fig. 4. Correlations of the remaining parameters with J[ and i for a binary with AM C Vn's parameter values as a function of cos i. 
Most correlations are weak for all inclinations except c$u. 
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Fig. 3. Two-dimensional error ellipses and cost extracted 
from the variance-covariance matrix, C, for HM Cnc. See 
Figure [2] for details. The 1-cr EM constraint of 7° in incli- 
nation reduces the corresponding 1-cr uncertainty in J{ from 
~ 1.22 X 1(T 23 to ~ 0.89 x KT 23 , roughly a factor of 1.4. 
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3.2. Influence of the inclination, i 

We have seen above that cos , depends on the inclination of the 
binary. For the special case of exactly face-on binaries (( = 0, 7r), 
the Fisher matrix is ill-defined owing to the symmetry in the sig- 
nal, namely, h(t, i + di) = h(t, t - di). The derivatives are zero, 
so the terms including l in the Fisher matrix are also zero, result- 
ing in an ill-defined matrix that cannot be inverted. In general, it 
is hard to track the behaviour of the variance-covariance matrix 
from the FIM alone. Hence, to determine the influence of incli- 
nation, we fix all parameters except M, i, and A, and calculate 
the matrices F and C for AM CVn as a function of t. This can be 
justified, as most correlations with J[ and t are not very strong; 
the largest correlation (of ~ 0.6) is for A as shown in Figure [4] 
This means that for all inclinations 



Cj\ cosi A 
Cfoffwnp 



Fig. 5. Signal-to-noise ratio, o- m , cr cosl , cr A , c Mc0Sl , c mA , and 
Ccosu as a function of cos t, keeping the other parameters fixed 
to the values of AM CVn. The vertical dashed line is the mea- 
sured value of t for this binary, 43°. The best known value for the 
amplitude is 9{ = 1.49 x 10~ 22 . As expected, the S/N is higher 
for face-on than edge-on orientations. The relative uncertainties 
in the amplitude and inclination are very large for inclinations 
up to 10°. The bottom-right panel shows the normalised correla- 
tions cj[ A as the dotted line and c cosl/1 as the solid line. 

thus we can focus on M, cos t, and A independently and reduce 
the 7x7 matrix to a 3 x 3 matrix, which is easier to interpret 
In Figure [5] we show the S/N, (relative) parameter uncertainties 

6 Another reason to consider a smaller matrix is that it is difficult to 
stabilise C for low inclinations i < 10°. We find that the optimal choice 
of d9j for which the parameter uncertainties tr, are stable varies as a 
function of inclination. For low values of i, the optimal d6j is orders of 
magnitude larger than for high values of i. Note that the FIM (as op- 
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Fig. 6. Signal-to-noise ratio and (relative) parameter uncertainties for AM CVn as a function of cos t. The general behaviour for 
cr^, <t C0S1 and cr A is consistent with Figure|5]above. The uncertainties in cr^, cr A , and crj behave as may be expected from the S/N, 
where parameters are more accurately determined for higher S/N. 



and the normalised correlation for these three parameters as a 
function of inclination. It is obvious that for a source at a given 
distance, the S/N for an edge-on binary is lower than that for 
a face-on binary owing to the less favourable orientation of the 
latter. Thus, the S/N is a function of cos t and has a minimum at 
l = 90° (cos t = 0). In general, we expect the parameter uncer- 
tainties to vary accordingly, i.e., a higher S/N should correspond 
to a higher precision and vice-versa. The longitude follows this 
behaviour in the bottom-left panel. However, for J[ and cos l this 
is not the case as the detector cannot accurately determine the 
inclination in the range of 0.7 < | cos i\ < 1, as shown in the top- 
right and middle-left panels. For i < 10° and l > 170° in partic- 
ular, the uncertainties in t and J[ become very large. However, 
as shown in the middle-right panel, the correlations are maximal 
for these almost face-on binaries, thus by constraining one of 
the parameters from EM (or other) data, we constrain the range 
of values that the other parameter can take. This method cannot 
be applied to edge-on binaries as the correlation is very weak in 
these cases. If we take a value of di for which the computation 
of the 3x3 variance-covariance matrix is numerically stable for 
all t and use it to produce the full 7x7 matrix, the parameter 
uncertainties behave as shown in Figure [6] The uncertainties cr$ a 
and cr^j behave similarly to cr^ and cr l since </>o, and if/ are also 
degenerate for nearly face-on systems. In contrast, the uncertain- 
ties erf, cr A , and <r sinj g behave as may be expected from the S/N 
since they are not strongly tied to any of the four parameters with 
strong correlations. 

3.3. The correlation between </> and if/ 

The orientation parameters <po and if/ are also highly anti- 
correlated with c M m -0.99 for AM CVn and HM Cnc. This 
can also be explained by examining the geometry of the phys- 
ical binary. The orbital phase <p is measured in the plane of the 
binary (perpendicular to the binary's angular momentum vector, 
L), while if/ is measured in the plane perpendicular to the line of 
sight, h. Both AM CVn and HM Cnc have an inclination in the 
range where the system appears to be face-on to the detector (see 
Section 3.1). Thus, </>o and if/ are practically degenerate for these 



posed to C) is stable for the same choice of dOj. In addition, the FIM 
choice of d6, that is made for the inclination affects only matrix ele- 
ments that contain terms related to i. However, this influence spreads to 
all matrix elements of C(= r~'), because of the matrix inversion. 



binaries. In contrast, for J0651 L and h are almost perpendicu- 
lar, and the effects of </>o and if/ can be more easily distinguished 
by eLISA. This correlation is not very noteworthy, because the 
parameters involved are of little astrophysical interest. 



3.4. The correlation between f and 4>o 

For J0651, there is a strong anti-correlation between / and <f>o, 
with a value of -0.89. This is because if the frequency of a sig- 
nal is changed by a small fraction, 6f, the signal accumulates a 
phase shift towards the end of the wave. This signal can be de- 
shifted in phase to obtain another signal with a similar match to 
the original, in turn making the two slightly dissimilar signals in- 
distinguishable. Hence, the two parameters / and </>o are degen- 
erate. This degeneracy is not present for AM CVn and HM Cnc 
owing to the strong anti-correlation between </>o and if/ and large 
errors in these two parameters. 



3.5. The correlation between f and f 

If we model HM Cnc with an additional eighth parameter, /, we 
find a strong anti-correlation between / and /. This can be ex- 
plained in a very similar way to the case of c/0 o above. Consider 
a signal which is obtained by changing the chirp by 8f, within its 
uncertainty. This will slightly change the frequency towards the 
end of the signal. By additionally changing / within its uncer- 
tainty, we can obtain a very similar signal to the one for which 
we changed only its chirp by 6f. Thus, the two parameters are 
degenerate for a mildly chirping binary. This could be a very 
useful correlation for constraining the / of the binary by using 
an EM constraint on /. However, the GW data analysis already 
provides a relative accuracy of ~ 10~ 6 for the frequency, which 
is hard to improve with EM observations. 

After providing a potential measurement of / for the J065 1 
source (M. Kilic, private communication), the correlation be- 
tween / and / could be used to test one of the predictions of 
General Relativity in Eq.|4] which holds for a binary that evolves 
only because of GW radiation loss. Given that the individual 
masses and their / are measured, Eq.|4]can be used to test the 
prediction of /. 
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4. Discussion 

The S/N's of the data for the verification binaries listed in Table[T] 
are relatively low. In particular, J0651 has S/N ~ 10. This may 
not satisfy S/N » 1, one of the criteria under which FIM studies 
are valid. To exclude the effects of S/N, we place the verification 
binaries at a closer distance so that they all have S/N of ~ 100. 
We find that all of the parameter uncertainties decrease as ex- 
pected, except for the uncertainty in M. However, the strong cor- 
relations we found for cj\ L , c^, Cf^ and c^y listed above in the 
matrices remain the same. Furthermore, these correlations do not 
change when we consider the six-link configuration of the clas- 
sic LISA with a five million km baseline and slightly different 
instrumental noise. From the six-link interferometer, three inde- 
pendent data streams can be synthesised, of which A, E, and T is 
one of the combinations. Naturally, this means that the S/N of the 
source will be higher, and that there may be additional informa- 
tion from these data streams. However, we found that when we 
consider these optimal TDIs, the strong correlations we found 
are unaffected. Additionally, the weak correlations between the 
rest of the parameters are also unaffected. 

It is well-known that in the low-frequency regime there are 
many overlapping sources, which can significantly degrade the 
parameter uncertainty. Crowder & Cornish (2004) performed a 
thorough FIM analysis and show that the degradation of param- 
eter uncertainties depends on the frequency, the number of bi- 
naries per frequency bin, and the observation time. However, the 
degradation is uniform for all the seven parameters, which means 
that each parameter is affected in the same way by the overlap- 
ping sources. Hence, the presence of overlapping sources affects 
the uncertainties in the parameters, but not the normalised cor- 
relation between the parameters. In addition the potential biases 
from imperfect subtraction of resolvable sources (Section 2.3) 
can affect the results, although this mostly affects the S/N degra- 
dation of the source and thus not the normalised correlations. 

Our study should be used to implement EM priors in an 
MCMC data analysis based on a Bayesian framework. While 
the latter leads to more accurate results, the FIM tools provides 
useful predictions ( VallisnerT]|2008| > . In the case of Galactic bi- 
naries, it has been shown that the one-dimensional marginalised 
posterior PDFs are well-matched by the FIM predictions when 
the data analysis of a single source is considered ( |Cornish "&] 
|Crowder|2005l . 

In a forthcoming paper, we will discuss the dependence of 
the aforementioned correlations and accuracies on the sky posi- 
tion and polarisation angle. Our test calculations show that the 
strong correlations c mi , c^, c/ A 
vary as a function of these parameters 

Finally, we observe that the uncertainty in the inclination of a 
binary indicates whether the system is eclipsing. In Section 3.1, 
we found that for J0651, the uncertainties in i, </>o, and tfr are small 
compared to those of HM Cnc, even though its S/N is a factor 
of three lower than that of HM Cnc. An uncertainty of cr cosi = 
0.043 translates to cr, ~ 2.5°. Hence, for J065 1 -like systems with 
S/N's of > 10, the small uncertainty inclination inferred from a 
GW analysis may be used to find candidate eclipsers. 



5. Conclusion 

We have performed Fisher-matrix studies to investigate whether 
there are correlations between the parameters that characterise 
Galactic binaries with short periods, between about six minutes 
and a few hours, which lie in the eLISA frequency band, 10~ 4 
Hz -1 Hz. We focused on three verification binaries, AM CVn, 



and Cif found above do not 



SDSS J065 1+2844, and the mildly chirping HM Cnc. Our main 
findings are: 

1. There are strong correlations between the strain amplitude 
and inclination, and between the phase and polarisation an- 
gle for AM CVn and HM Cnc. For the latter, there is an ad- 
ditional strong anti-correlation between / and / when it is 
modelled as a mildly chirping source. 

2. The remaining parameters are not or only weakly correlated 
for the binaries considered. 

3. The correlation between strain amplitude and inclination can 
be very useful in constraining the amplitude of the binary, 
which is a function of the two masses, the distance to the 
binary, and its frequency. Since the frequency is very accu- 
rately determined (of the order of \0 9 Hz) by the GW data 
analysis, a combination of the masses and the distance can 
be constrained by using an EM constraint on the inclination. 

4. These correlations depend strongly on the inclination of the 
system: approximately face-on (t < 45°, l > 135°) binaries 
have strongly correlated parameters, whereas for edge-on bi- 
naries, the correlations become very weak. 

5. For binaries with very low inclinations (i < 10° and t > 
170°), the uncertainties in amplitude, inclination, phase, and 
polarisation angle derived from GW data are very large. 

6. We have found that the influence of S/N on the correlations 
is not significant; placing the sources at shorter distances to 
give them higher S/N yields the same normalised correla- 
tions. 

7. The strong correlation between amplitude and inclination is 
the most useful correlation found, since inclinations are typi- 
cally measured to a higher accuracy by EM observations than 
with GW measurements alone. Hence, this correlation can be 
used to constrain the amplitude by a factor of six for systems 
with orientations and S/N's similar to those of AM CVn. 

8. We have found that some correlations also depend on sky 
position and polarisation angle, which we will discuss in a 
forthcoming paper. 

9. Even for signals with an S/N of ~ 10, we have been able to 
reliably determine the inclination when the system is edge- 
on. This will enable efficient searches for eclipsing systems 
with EM instruments. 
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Appendix A: Stabilisation of the 
variance-covariance matrix 

Here we give an example of the Stabilisation of the variance- 
covariance matrix for AM CVn. The dQ; are varied for each 
parameter over seven orders of magnitude as described in 
Section 2.4. In this example, d9j\, dOf, d0 L , dd^, are varied from 
(10~ 4 , 10 2 )x yfexOj, where, 0,- are the true values for AM CVn. 
The values of dO^, d9p, d9,\ are varied from 10~ 6 , ... , 10°, more 
than for parameters above, because the numerical derivatives of 
<p, ft, and A are unstable for lower values . The first row of the un- 
certainties in the parameters correspond to the first set of d6u et 
cetera. In principle, the uncertainties cr,- should level off around 
the third row. The uncertainties for all our verification binaries 
are calculated in this way. 




f[Hz\ 

Fig. B.l. Average PSDs of instrumental noise for eLISA with 
varying sampling time, At. For lower-frequency sources such 
as AM CVn and J0651, a sampling as fine as 64s can be used. 
However, for HM Cnc, whose frequency is higher, a high At 
would overestimate its S/N, since the level of the noise is too 
low for the longest sampling times. 



Table A.l. Parameter uncertainties, <x,, in a range of d8, for the 
case of AM CVn, with S/N ~ 1 1 for T \, s = 2 years. The values 
start to level off from the third to fifth row. 
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Appendix B: Sampling 

In Fig |B.1| we show the effect of the sampling time At on the 
noise PSD of the instrumental noise, which is the average spec- 
trum of the TDI X observable. At / < 5 x 10~ 3 Hz, the sampling 
time does not affect the level of the PSD. Above these frequen- 
cies, the PSD with a relatively long sampling time of 64s under- 
estimates its true level. 
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